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Abstract 

We determine the leading-order nonrelativistic quantum chromodynamics (NRQCD) matrix el- 
ement (Ci)t and the ratio (q 2 )-r, for T = T(nS) with n = 1, 2, and 3 by comparing the measured 
values for T[T — > e + e~] with the NRQCD factorization formula in which relativistic corrections 
are resummed to all orders in the heavy-quark velocity v. The values for (q 2 )r ; which is the ratio 
of order-?; 2 matrix element to (Oi)x, are new. They can be used for NRQCD predictions involving 
T(nS) and r] b (nS) with relativistic corrections. As an application, we predict the two-photon decay 
rates for the spin-singlet states: T[rj b (lS) -> 77] = 0.512+°;°^ keV, T[r] b (2S) 77] = 0.235+^;^ 
keV, and T[r] b (3S) 77] = 0.170l£^ keV. 



I. INTRODUCTION 



The pseudoscalar bottomonium rjb(lS), which is the spin-singlet S'-wave ground state, was 
first observed in the photon energy spectrum of the radiative T(3S) decay [1] and confirmed 
in the radiative T(2S) decay [2] by the BABAR Collaboration. The state was also confirmed 
by the CLEO Collaboration again in Y(3S) — > •yrjb^lS) [3]. So far, only the mass for the 
rjb(lS) is known as m Vb (i S ) = 9390.9 ± 2.8 MeV [4], and any of its exclusive decay modes 
has not been observed, yet. Among its various decay modes, recent theoretical studies have 
been concentrated on relatively clean channels like r/ b — > J/ip J/t/j [5-8], r] b — > J/ip 7 [9], and 
others [10, ll]. 1 On the other hand, the most elementary exclusive decay channel is r] b — > 77, 
although it has a large background. With the decay mode T — > e + e~ of the spin-triplet 
partner, r) b — > 77 must be well described by the nonrelativistic quantum chromodynamics 
(NRQCD) factorization formulas for the electromagnetic decay of heavy quarkonia [12]. If 
one makes use of the heavy-quark spin symmetry, then one can make a rough estimate of 
the decay rate, whose branching fraction is ~ 10~ 5 , which is relatively greater than other 
channels listed above. 

Available predictions for the decay rate T[r] b — > 77] are based on the potential model [13- 
19], the Salpeter method [20-22], or the heavy-quark spin symmetry [23]. Some of them in- 
clude the effects of the relativistic corrections and binding effects and most of the predictions 
rely on the heavy-quark spin symmetry between the spin-singlet and -triplet states. One 
can estimate the spin dependence of the rate systematically by making use of the potential 
NRQCD [24, 25]: In Ref. [26], the decay rate was computed to the next-to-next-to- leading 
logarithmic accuracies as T[r] b {lS) 77] = 0.659 ± 0.089 (th.)+°;^ (5a s ) ± 0.015 (exp.) 
keV. Recently, an updated potential-NRQCD prediction for the decay rate became avail- 
able: 0.54 ± 0.15 keV [27], in which leading relativistic corrections are included. 

In the mean time, there has been a significant progress in the NRQCD calculations for S'- 
wave charmonium production and decay, in which relativistic corrections of all orders in the 
heavy-quark velocity v are resummed [28] . Precise determination of the wavefunction at the 
origin for the J/t/j was made based on this method [29-32]. This method has been applied 
to reconcile the large discrepancy between the theoretical prediction and the experimental 

1 Throughout this Letter we use collective notations T]b and T that indicate rib(nS) and T(nS), respectively, 
for n = 1, 2, and 3 unless a specific state is specified. 
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results for the cross section a[e + e~ — > J/ip + i] c ] at the B factories [33-35]. Therefore, it 
is worthwhile to improve the NRQCD prediction for T[r]b — > 77] by taking into account 
the relativistic corrections to all orders in v. In order to carry out such an analysis, one 
needs to know the values for the color-singlet NRQCD matrix elements and those involving 
relativistic corrections. The NRQCD matrix element for the S-wave bottomonia can be 
determined by making use of the measured values for r[T — > e + e~] up to corrections of 
spin-symmetry breaking effects. Unfortunately, available order-f 2 NRQCD matrix element 
that has been fixed from lattice QCD simulations [36, 37] suffers from large uncertainties 
originated from slow convergence of the cut-off regularization method. 

In this Letter, we first determine the NRQCD matrix elements for the .S'-wave bottomo- 
nium states that are required to compute the relativistic corrections with considerably less 
uncertainties than available values, extending the method in Refs. [28, 30]. As an application, 
we compute T[rj b — > 77], in which corrections of order the strong coupling a s and relativistic 
corrections of all orders in a° s q 2n are included. Here, q is half the relative momentum of b 
and b in the bottomonium rest frame. The remainder of this Letter is organized as follows: 
In Section II, we estimate the NRQCD matrix elements for the S'-wave bottomonium states 
by making use of the resummed NRQCD factorization formula against empirical data for 
the spin-triplet states. Our prediction for r[rii> — > 77] is presented in Section III with the 
comparison with available predictions and we summarize in Section IV. 

II. NRQCD MATRIX ELEMENTS FOR THE T(nS) 

In this section, we briefly review the method to determine the NRQCD matrix elements 
for T at leading and subleading order in q 2 based on the strategy for the charmonium 
counterpart in Ref. [30]. The results are compared with those of lattice QCD calculations. 

The NRQCD factorization formula for the electromagnetic decay of the S'-wave quarko- 
nium if is a linear combination of nonperturbative NRQCD matrix elements {O n )u that are 
classified in powers of v, where O n is the NRQCD operator. The factorization is achieved 
at the amplitude level and the ratio (q 2n )H of the order-qr 2n matrix element to the leading 
one are all, in general, independent. In addition, the ratios (q 2n )H have power-ultraviolet 
divergences that must be regulated and, therefore, the values can even be negative under 
subtraction. In lattice QCD calculations, this subtraction is made by making use of the 
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hard-cut-off regularization whose convergence is slow, resulting in large uncertainties [28]. 
However, in an electromagnetic decay, in which the color-singlet contributions dominate, 
one can calculate the quarkonium wavefunction of the leading heavy-quark-antiquark (QQ) 
Fock state up to corrections of relative order v 2 if one knows the static, spin-independent 
QQ potential exactly. The authors of Refs. [28, 30] have constructed the generalized ver- 
sion, {q 2n }n = [(q 2 )H] n , of the Gremm-Kapustin relation [38] to resum a class of relativistic 
corrections. The method has been devised to be consistent with dimensional regularization 
of these power-ultraviolet divergent matrix elements. 

The resultant formula for the decay rate of T — y e + e~, in which relativistic corrections 
of all color-singlet QQ operator matrix elements are resummed, is given by [29, 30] 

Vi>r, (1) 

where mj is the T mass, (Oi)r is the color-singlet NRQCD matrix element for the electro- 
magnetic decay of the T at leading order in v, and (v 2 )x = (q 2 )r /m 2 with the bottom-quark 
mass m&. The resummed relativistic corrections to all orders in v at order a 2 a° s are con- 
tained in the function f(x) = x/[3(l + x + y/1 + x)} with x = (f 2 )r and in the factor l/m\ 
implicitly. 

The order-a 2 corrections to T[T — y e + e~] (Refs. [39, 40]) contain a strong dependence on 
the NRQCD factorization scale. If one were to include those corrections in Eq. (1) and use 
it to determine (Ci)x, then (0i)r would also contain a strong dependence on the NRQCD 
factorization scale, which would cancel in other quarkonium decay and production processes 
only if the short-distance coefficients were calculated through relative order a 2 . Generally, 
short-distance coefficients for quarkonium processes have not been calculated beyond relative 
order a s . For this reason, we omit the order-ct 2 corrections to the leptonic width in Eq. (1). 
Nevertheless, if one includes the order-a 2 corrections and take the factorization scale to be 
m b , the resultant NRQCD matrix elements are increased by about a factor of 40%. 

We briefly discuss the method employed in this Letter to compute (Oi)x and (q 2 )r- 
We follow the method given in Ref. [30] and make use of the Cornell potential model [41]. 
By using the Schrodinger equation we can express (Oi)r and (q 2 )r as functions of the 
parameters of the Cornell potential model, which are the mass parameter in the Schrodinger 
equation, the string tension, and the Coulomb strength of the Cornell potential. The mass 
parameter can be expressed in terms of the 1S-2S mass splitting [30], which we compute 
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from the masses of T(IS') and T(2S'). The value of the string tension, which is universal, 
is taken from lattice measurements as 0.1682 ± 0.0053 GeV 2 [42-45]. Finally, the Coulomb 
strength parameter is determined by constraining the rate (1) to be consistent with the 
experimental value [4] and solving the resulting nonlinear equation numerically. Because of 
this, the value of the Coulomb strength parameter is chosen differently for each quarkonium. 
From the fixed values of the model parameters we obtain the numerical values of the matrix 
elements. For details of the method, we refer the readers to Ref. [30] and references therein. 

We list the numerical values and uncertainties of the parameters used in Eq. (1). The 
measured leptonic widths of T(nS) are T[T(1S) -> e+e"] = 1.340 ± 0.018 keV, T[T(2S) ->■ 
e + e -] = 0.612 ± 0.011 keV, and r[T(3S) ->■ e + e~] = 0.443 ± 0.008 keV [4]. The masses 
for the T(nS) states are taken to be tot(is) = 9.46030 GeV, m T (2S) — 10.02326 GeV, and 
m x(35) = 10.3552 GeV [4], where the errors (< 5 x 10 _3 %) are neglected. The factorization 
formula (1) depends on m b implicitly through {v 2 )r, where we use the one-loop pole mass 
mj = 4.6 ±0.1 GeV. We evaluate a(/i) and a s {ji) at the scale, the momentum transfer at 
the quarkonium-photon vertex. The values are a(/x) = 1/131 in every case, a s [mr(is)] = 
0.180 ± 0.032, a s [m r( 2S)] = 0.177 ± 0.031, and a s [m r{3S )} = 0.176 ± 0.031, where the 
uncertainties of relative order a s are included in the strong coupling. The main difference 
between this analysis and that for the S'-wave charmonium in Ref. [30] is that there are no 
measured data for T[i] b — > 77]. Therefore, we use the spin-triplet data only. 

By carrying out these calculations, the Coulomb strength parameter is fixed as 9.955 
for T(1S), 10.960 for T(2S), and 11.127 for T(3S), respectively. From these we ob- 
tain our results for (Oi)r and (q 2 )r, which are tabulated in Table I. The correspond- 
ing values for the quantity (v 2 ) r are (v 2 ) r{is) = -0.009^;^, (v 2 ) r{2S ) = 0.090±g:oii. and 
(v 2 )r(3S) = 0.155^Q;o}g. These values are in rough agreements with the typical estimate 
v 2 ~ 0.1 for the bottomonium except that (v 2 )r(is) is tiny. The error bars in Table I reflect 
the uncertainties arising from m&, r[T — > e + e~], string tension, a s , and the ignorance of the 
spin-dependent interactions of the potential in the Schrodinger equation [30], all of which 
are added in quadrature. The values for the leading-order matrix elements (0±)r in Table I 
have been used to predict the inclusive charm production in T decays [46] and the heavy 
quarkonium production associated with a photon in e + e~ annihilation [47]. The values for 
(q 2 )r in Table I are new. The value for T(IS') has errors significantly less than those of 
the available lattice QCD calculations [36, 37]. One can reduce theoretical uncertainties 



5 



TABLE I: The NRQCD matrix element (Oi)x at the leading order in v in units of GeV 3 and ratios 
(q 2 ) r in units of GeV 2 for T = T(1S), T(2S), and T(3S). 



Sources of errors 


(Ci)t(is) 


(q 2 )t(is) 


(Cl)T(25) 


(<? 2 )t(2S) 


(Cl)T(35) 


(Q 2 )t(3S) 


Arrib 
others 


3.069^;^? 
3.069^ 


-n ioq+o.ooo 

U. iyO_Q QQQ 

n iqq+0.072 
u.iyo_ 073 


1 RO9+0.002 
l.OZO_ 002 

1 623+ - 112 
±.ozo_ 103 


1 saa+ - 001 
l.»y8_ .000 

1 SQS+ - 210 

i.8y»_ 0210 


-, 97Q +0.003 
i.z/y_ 003 

1 O7Q+0.090 
L - Zl y -0.083 


q ooq+0.003 
O.ZOO_ .002 

q ooq+0.353 

a.zoa_ .352 


total 


3.06918;™ 


u.±»o_ 073 


1 623+ - 112 
±.uzo_ 103 


1.898«; 2 i° 


i 27Q+ - 090 
i.z/y_ .083 


q ooq+0.353 

O.zoo_ 352 



by considering the dependence on m b that also appears in the short-distance coefficients of 
factorization formulas. Therefore, we present the sources of errors in Table I. Unlike the 
S-wave charmonium case in Ref. [30], the uncertainties of (Oi)r and (q 2 )r due to the errors 
of the heavy-quark mass are insignificant. 

Our results are now compared with those for the ground-state S'-wave bottomonium 
obtained from a lattice QCD simulation. The results from the quenched approximation are 
given in Ref. [36]. We quote the updated results of the unquenched analysis in Ref. [37]: 
The leading-order NRQCD matrix element is (Oi) ls = 4.10(1)(9)(41) GeV 3 and the ratio 
(q 2 )is ranges from about —5 GeV 2 to about 2 GeV 2 [37]. Here, the subscript IS indicates 
the average of 7)1,(1 S) and T(1S). Our central value for the (Oi)t(is) is about 25% smaller 
than that of Ref. [37], which is greater than the quenched case [36] by about a factor of 2. In 
the case of (q 2 )y(is), our result is consistent with that in Ref. [37] but ours has uncertainties 
significantly smaller than that of the lattice result. 

III. TWO-PHOTON WIDTHS FOR THE t) b 

In this section, we predict — > 77] by making use of the NRQCD matrix elements 
determined in Section II. In fact, the NRQCD matrix element (Oi) Vb that appears in the 
factorization formula for T[7) b — > 77] might be different from (Oi)r by a relative order v 2 , 
which breaks the approximate heavy-quark spin symmetry [12]. We recall that the effect of 
spin-symmetry breaking in the low- lying S'-wave charmonia J/ip and 7) c is not significant [30]. 
Therefore, the errors in the approximation (Oi) Vb ( n s) ~ (Ci)x(nS) ma y be insignificant based 
on the fact that (v 2 ) T < (v 2 )j/^ = 0.22 [30]. 

As in the leptonic decay of the T, we include the relativistic corrections to T[7) b — > 77] to 
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TABLE II: Two-photon widths of the rjb(nS) in units of keV. 



State 


T) b (lS) 


rj b (2S) 


%(35) 


r 

i 77 


n n n+0.096 

u.o±z_ 0094 


n 90C+0.043 

U.ZOO_ 043 


17 q+0.031 
u - i ' u -0.031 



all orders in v. The resultant factorization formula is given by [30, 31] 

2 

(Oi)*, (2) 

where the relativistic corrections are incorporated into the function g(x) = 1 — {log[l + 
2^/xil+x) + 2x]}/[2y/x(l + x)] with x = {v 2 ) Vb = {q 2 ) Vb /ml. The input parameters 
for the numerical calculations are chosen in a similar way in Ref. [30] for T[i] c — > 77]. 
The scale /i for the couplings a and a s are taken to be the momentum transfer at the 
photon-heavy-quark vertex, namely, m Vb /2: a = 1/132 for every case, a s [m Vb (is)/2] = 
0.216 ± 0.046, a s [m Vb{2S) /2] = 0.212 ± 0.045, and a s [m rjb{zs) /2\ = 0.210 ± 0.044, where 
the uncertainties of the strong coupling are of relative order a s . For the meson masses, 
we use m m[lS ) = 9390.9 ± 2.8 MeV [4], m Vb{2 s) = 9.97 GeV, and m m{?S ) = 10.3 GeV, 
where we have assumed that mx(„s) — m Vb ( n s) = 0.5 MeV for n = 2 and 3. While this 
value for the hyperfine mass splitting is smaller than the measured value for the IS states 
m x(i5) — m r] b (is) — 69.3 ± 2.8 MeV [4], it is comparable to that for the 2S charmonia, 
m ip(2S) — m i]c(2S) — 49 MeV [4]. Note that the uncertainties from m^^s) an d m %(3S) are 
insignificant because the factorization formula (2) does not depend on them but on m^. Like 
the leptonic width of the T [Eq. (1)], we omit the order-a^ corrections to the two-photon 
width of the rjb, whose result is available in Ref. [48]. 

The resultant predictions for — > 77] are tabulated in Table II. The errors include 
the uncertainties of a s , m b , and the values for (Ci)t and (q 2 )x in Table I. We also include 
the errors of using (Oi)y and (q 2 )r, which are of relative order v 2 set to be 0.1. From 
the order-a 2 corrections to the electromagnetic widths of the T and r/ b [39, 40, 48], we find 
that the order-a 2 corrections account for —2.64a 2 in the ratio Tf^ — >■ 77]/T[T — > e + e~], 
if we choose the NRQCD factorization scale to be mb- Therefore, we include the errors of 
omitting the order-a 2 corrections in using T[T — > e + e~] to determine Tfr^ — > 77] as 2.64a 2 . 
This implies that the large correction to the leading-order NRQCD matrix elements arising 
from inclusion of the order-a 2 corrections, as briefly shown in the previous section, almost 
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cancels the order-a 2 corrections to the two-photon width of the rj b . All of the errors listed 
above are added in quadrature. We can compare our results with previous predictions. In 
the case of 7] b (lS), available predictions range from 0.170 keV to 0.659 keV. The results in 
Refs. [15-18, 20-23] agree with our prediction within errors, while some models [13, 14, 19], 
which does not use the heavy-quark spin symmetry, apparently underestimate the rate in 
comparison with ours. Our result T[r]b(lS) — > 77] = 0. 512^0^94 agrees with the most recent 
potential-NRQCD prediction 0.54 ± 0.15 keV in Ref. [27] in which the leading relativistic 
corrections are included, while it is smaller than another potential-NRQCD prediction in 
Ref. [26]. Note that we have borrowed (Ci) T and (q 2 )x for (Oi) Vb and (q 2 ) Vb in Eq. (2) 
after taking into account the errors of spin-symmetry breaking effect as v 2 ~ 0.1 because 
r[?7ft — > 77] are not measured. Once T[rj b (nS) — > 77] are measured in the future, one can 
determine (Oi) Vb ( n s) and (q 2 ) Vb (nS) (° r eventually (v 2 ) Vb ( n s)) with an improved accuracy in 
combination with the measured values for T[T — > e + e~]. 

IV. SUMMARY 

In summary, we have determined the leading-order NRQCD matrix element (Ci)x and 
the ratio (q 2 )r, for T = T(nS) with n = 1, 2, and 3 by comparing the measured values 
for the leptonic decay rates of the T with the NRQCD factorization formula in which rel- 
ativistic corrections to all orders in v are included. The values for (q 2 )r are new and can 
be used for various phenomenological predictions for T and rj b including relativistic cor- 
rections. The values for (q 2 )y are consistent with the naive expectation of the velocity- 
scaling rules except that (q 2 )r(is) is tiny By assuming approximate heavy-quark spin 
symmetry with the uncertainties of relative order v 2 ~ 0.1, we used (Oi)x and (q 2 )r 
to estimate T[r] b (lS) 77] = 0.512^°;°^ keV, T[r] b (2S) 77] = 0.235^;^ keV, and 
T[rj b (3S) 77] = 0.170^03} keV. Our prediction for T[r} b (lS) -> 77] is consistent with a 
recent potential-NRQCD prediction in Ref. [27], in which the leading relativistic corrections 
are included. 

By making use of the ratio T[rj b (lS) — > r yy]/T[ri b (lS) — > gg], one can make a rough esti- 
mate of the branching fraction for t] b (lS) — > 77 as ~ 6.9 x 10^ 5 . The BABAR Collaboration 
reported 19200 ± 2000 ± 2100 rj b (lS) events out of (109 ± 1) x 10 6 T(3S) samples [1]. They 
also obtained 12800 ± 3500±^fgg r] b (lS) events from (91.6 ± 0.9) x 10 6 T(2S) samples [2]. 
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These are not sufficient to observe the mode r] b (lS) — > 77. However, we expect that this 
channel can be observed at the superKEKB or superB factory if more data are accumu- 
lated. The CERN Large Hadron Collider is expected to produce about 5 x 10 9 rj^s with 
the integrated luminosity ~ 300 fb _1 [6], with which one can probe about 45000 events of 
rjb(lS) — > 77. We anticipate such a stage against which our predictions can be tested. 
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